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Abstract 

We present the complete solution of the first order metric and 
density perturbation equations in a spatially flat {K = 0), Friedmann- 
Robertson- Walker (FRW) universe filled with pressureless ideal fluid, 
in the presence of cosmological constant. We use covariant linear 
perturbation formalism and the comoving gauge condition to obtain 
the field and conservation equations. The solution contains all modes 
of the perturbations, i.e. scalar, vector and tensor modes, and we 
show that our results are in agreement with the Sachs & Wolfe metric 
perturbation formalism. 
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1 Introduction 

Recent cosmological measurements by the Wilkinson Microvawe Anisotropy 
Probe (WMAP) P], la type supernovae observations [2j and other sky sur- 
veys eg. Sloan Digital Sky Survey (SDSS) [H] seem to prove today that the 
expansion of the universe is accelerating. Although there exist many differ- 
ent approaches explaining this phenomenon, the most accepted description 
in the literature is the concordance model. In this model the universe is fiat 
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{K = 0), homogeneous and isotropic with nonzero cosmological constant A 
in the Einstein equations. 

In this work our aim is to obtain the complete first order solution of 
the perturbations of a concordance model filled with pressureless matter in 
Bardeen's covariant formalism Hu & Sugiyama j3] claim that in the pres- 
ence of cosmological constant one must use a numerical approach to get the 
result. Recent calculations of Perjes et al. [U] yield all the C°° solution of 
the problem following the original work of Sachs & Wolfe |2|- Here we show 
that the complete analytic solution of the perturbed field equations can also 
be obtained using the covariant linear perturbation formalism, in agreement 
with the results of Perjes et al. jSj. Providing the solution of the problem 
in this formalism has the advantage that recent model calculations of the 
fluctuations of the Cosmic Microwave Background Radiation (CMBR) pre- 
fer using Bardeen's covariant linear formalism rather than the perturbation 
approach of Sachs & Wolfe. 

In Sec. 2, we present the well known solution of the background quantities 
in a fiat, homogeneous and isotropic universe and give a brief description of 
the covariant linear perturbation formalism. The perturbed field equations 
and their solution for scalar, vector and tensor modes are presented in Sec. 
3. In Sec. 4 we build up the complete first order solution of the metric and 
show the agreement with the results of [B]. 

2 Background metric and linear perturbations 
in covariant formalism 

As background solution we use the homogeneous and isotropic FRW metric 



where 7^j, is the three metric of a space with constant spatial curvature 
K, and rj is the conformal time variable. Introducing the comoving time 
coordinate t by the relation adi] = dt, the scale factor a(t) can be expressed 
in the case of a flat {K = 0) space and pressureless {p = 0) matter source 
as IE] 
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Throughout this paper we use units in which the gravitational constant G = 
I/Stt, the speed of hght c = 1 and we set Cq = 0. Roman indices run from 
to 3 and Greek indices run from 1 to 3. We use 7^,^ and its inverse to raise 
and lower spatial indices of first order quantities. 

To obtain the perturbed field equations we use Bardeen's covariant linear 
formalism |3] as it is presented by Hu [H] • Hereafter we discuss only the main 
steps of this description to derive the basic equations. 

For a spatially flat FRW universe the general perturbations of the metric 
tensor have the following form 

^00 = -a\l + 2A) , 

goa = —a^Ba , (3) 

The functions A, Hi, and HTap give a complete representation of the 
metric, where Hxap is a 3 x 3 trace- free tensor. 

We define an effective energy momentum tensor in the Einstein equations 
with cosmological constant as 

The energy momentum tensor Tik in the case of a homogeneous and isotropic 
ideal fluid source can be generally perturbed as 

-p-6p , 

{p + p){va-B^) , (5) 

-(p + pK , 

where 6p, 6p and Va are the density, pressure and velocity perturbations, 
respectively, and we introduce the effective energy density p and pressure p 
as new variables with the following definition 

p = p + A, p = p-A. (6) 

The n"g tensor represents the anisotropic stress perturbations. 

In covariant linear formalism the Einstein equations can be decoupled 
into a set of ordinary differential equations by employing scalar, vector and 
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tensor eigenmodcs of the Laplacian operator which form a complete set. In 
a spatially flat {K — 0) universe, these eigenmodes are plane waves 

Q(°) = exp(ik-x), 

Qi^'^ = ^(ei±^e2)aexp(zk-x), (7) 



/ 3 

Q'op^ = -J-{ei±te2)a{ei±ie2)pexpiik-x) , 



with the unit vectors ei and e2 spanning the plane transverse to the wave 
vector k. For an arbitrary scalar, vector and tensor function the components 
of the kth eigenmode become 

F(x) = F(fc)g(°), F„(x)=X: F^"^\k)Q^r\ ^a/3(x)=E F^"^\k)Q%\ 

m=—l m=—2 

(8) 

In the plane wave expansion there are relations between the curl free vectors 
and longitudinal components of tensors, and the scalar and vector modes of 
the eigenfunctions as follows 

= (A:-V„V, + ^7a/3)g(°\ (9) 

Q?/^ - -^[v.Qr + v,g(-)], 

where V is the covariant derivative operator with respect to ja/s- 

In the perturbed space-time we choose comoving coordinates. Using this 
gauge the following conditions hold 

A = , t;" = . (10) 

In further calculations we are interested in a spatially flat (K = 0) universe 
filled with pressureless {p = 0) ideal fluid and vanishing anisotropic stress 
perturbations (11"^ = 0). 



3 The field equations and their solution 

In this section we solve the Einstein equations for scalar, vector and tensor 
modes. First we present the field equations using the conditions discussed at 
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the end of the previous section, then we give the complete solution for each 
mode of the perturbations. 



Scalar modes 



The field equations for scalar modes have the following form 



, B 

Hl + -Ht + a[ — 



1 



a 
~2 



5p — 3dp— 

k . 



d 



rr 1 z> ^ 

aHL + -aHT = —p— 



d : 

a— + a 
dt 



) {kB 



aHr) = 



(11) 



' z> kB' 
aHL + — 



6 ^ 



and the scalar modes of the conservation equations become 



d 



d 



a— + 3d \ Sp = —SpaHi , a— + 4a p 



dt 



dt 



B 



. 



(12) 



For the scalar quantities the solution of the field equations gives 
Sp 



'""'^^^'^ [K,ik) - K,{k)m] , B ^ 



sinh^(Ct) 

^3 



Cm cl 



^ coth(Ct) [K^{k) - K2{k)I{t)] + V{k) 



M 



-3^i-^^coth(Ct) + i/o(A:) , 



(13) 



where Ki{k), K2{k), Bo^k), V{k) and HQ{k) are arbitrary scalar functions 
depending on spatial coordinates. The function I{t) is the following integral 
of the time variable 



I(t) = 2-2/=^v^ / 

JO 



sinh2/^(C 



cosh (Cr) 



dr 



(14) 



The Legendre normal form of the elliptic integral I{t) is given in the appendix 
of 0. 
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Vector modes 



The Einstein equations for vector modes are 



5(±i) 



a— + za 

at 



- aH, 



(±1) 



(15) 
(16) 



and the conservation equations become 



a— + Aa 
at 



5(±i) 



(17) 



Equations (fTK|) - (fTTjl are not independent. Inserting Eq. (|TH|l into Eq. (fT?)|) we 
get Eq ifTTj) . The general solution of this system of equations is 



5(±i) 



k R 



(±1) 



Cm 0, 



(18) 



H, 



(±1) 



5, 



(±1) 



^2 2 

— J(t) coth(Ct) 

Cm ao 



(±1) 



where BQ^^\k), and HQ^^\k) are arbitrary functions depending on only spa- 
tial coordinates. The function J{t) is also an elliptic integral of the time 
variable t, and can be expressed in terms of I{t) as follows 



■m 



^/(t) -3sinh-^/3(Ct)cosh-^(Ct) 



21/3 



(19) 



Tensor modes 



For tensor modes we get a source-free gravitational wave propagation equa- 
tion 

a2^f(±2) + 3aaH^^'^ + k'^H^^''^ = (20) 

as is the case in the absence of cosmo logical constant. The general solution 
of the wave equation in the presence of A is given in Eq.(50) of jH]. 
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4 The complete form of the metric 

Having in hand the solution for all modes of the perturbations, Eqs. fll3p 
and ()18|) , we can build up the general form of the metric (jSj) . The first order 
part of the metric takes the following form 



(+2)^(+2) 




(-2)^(-2) 

:/3 



2 

ao 



- — coth(Ct) + H^^^^ 



(+1) 



— coth(Ct) + 



(-1) 



ao 



Wo 



+ 



Cm 
Bo 



^ coth{Ct) [Ki - K^Iit)] -W - — ^ coth(Ct) + 



2Cal 



Ho 



2Cal 3 



coth(Ct) - 



7a/3 \ Q 



(0) 



and 



ka 



9oo 



C, 



M 



BoQf+B^^'^Q^^'^+B^-'^Q 



(-1) 



kakfS 

~k^ 

(21) 
(22) 



To show the equivalence of the metric (PT|) and (f^^ with the C°° solution in 
[H], we choose the space dependent integrational functions as follows, 



Bo 
Ho 



ec 



M 



2al 
—kCMCf^ 



3eC 



M 



27/3a5C'2 



B, 



0, B'o 



(±1) 



-kCMCl"^^^ 



(23) 



Ht'^ = . 



Here Ak, Bk, Cl and C^ are the components of the fcth eigenmode of 
the spatial functions in Eqs. (52-54) of [H]. Inserting the values in Eq. ()23j] 
into Eqs. (PT|) and we find the agreement of the two solutions up to the 
remaining gauge freedom, which is not completely fixed in [H]. 



5 Concluding remarks 

The result of our work is relevant in the context of studying the CMBR fluc- 
tuations in the picture of first order approximation of a concordance model. 
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We presented the complete form of the metric perturbations with cosmolog- 
ical constant in Bardeen's covariant hnear formahsm. We concluded that all 
the analytic solution can be obtained in agreement with the Sachs & Wolfe 
approach. This result can be useful in predicting the power spectrum of the 
fluctuations in the microwave background radiation. 

We dedicate this paper to the memory of our supervisor Prof. Zoltan 
Perjes who took part actively in the first period of the calculations. Regret- 
tably we had to finish writing the manuscript without his help and sugges- 
tions. 
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